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Why the interest in fusion propulsion? 
Go Far @ Get There Fast @ Carry Lots of Stuff 


Steady state 


Quiet clectric discharge 12 5,500 500, — — 83 

IEC 0.02 3,000 300 ae 8a UF eve 

Gas dynamic mirror 133 142,000 1,500 2.2 x 10% —_—— 85 Mo 

Tandem mirror 2.2 =- = 450 5.0x10” —- 86. 

Spheromak 5.75 20,000. 1,033 8.0 x 102 — 87 

FRC 1 105-108 1,100 1.0.x 10?! —— 68 * Solar system destinations 

Colliding beam FRC S. 1.4 x 10° 33 5.0x 10” —_—— 88 F 4 5 

Dipole 1 104 1.300 1.0x 102 ——_. 89 require Av = 10*— 10° m/s 

Spherical torus 8.7 35,435 1,690 5.0x 10” _—— 66 eu a ; 

Spheromak 0.375 — 10,000-40,000 540 (dry) 4.0.x 107 —_—— 67 High lop ( v./8) required 
Pulsed for reasonable payload 

Inertial fusion rocket | 270,000 (ays 1.0x 10% 100 87 mass fractions 

Inertial confinement fusion 6.62 27,200 6,000 1.0 x 10? 30 76 

MTF 16.8 70,485 121 1.0.x 107 20 48 * Fusion offers high specific 

Pulsed high density FRC 16 50,000. 67 2.0 x 102 10 69 power (kW/kg) to reduce 


Pulsed Z pinch 16.2 19,436 598 1.0.x 1078 10 82 SHE 
trip times 


Cassibry et al., “Case and Development Path for Fusion Propulsion,” Journal of Spacecraft and 
Rockets, 52 (2), pp. 595-612, March-April 2015 (and references therein) 


Several studies have defined the benefits of fusion propulsion for rapid trip times 
to Mars and the outer solar system 
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Magnetoinertial Fusion (MIF) 


¢ Ahigh current pulse generates an intense magnetic field that 
implodes an annular liner, heating and confining a plasma on axis 


Lindl-Widner Diagram (DT cylinder at stagnation)* 
100 


¢ A strong internal magnetic field 
reduces electron thermal losses 
and enhances alpha particle 
self-heating of the fusion fuel 


* Reduces the areal density (pr) 
threshold for ignition to well : 
below conventional Inertial : at 
Confinement Fusion values 


fuel temperature (keV) 


: . 
* Currently used in liner-driven — " 


fuel pR (g/cm’) 


FRC and Z-pinch experiments 


*Basko, M. et al., “Ignition conditions for magnetized target fusion 
in cylindrical geometry,” Nuclear Fusion, Vol. 40, No. 1 (2000) 


Liner Dynamics 


Paths to MIF Compression 


¢ Z-pinch: Pulsed axial current through liner generates azimuthal 
magnetic field; radial j,.XB, Lorentz force implodes the liner 


* ©-pinch: Induced azimuthal liner current generates axial 
magnetic field; radial j.XB, Lorentz force implodes the liner 
¢ Liner material serves to: 
- Conserve magnetic flux 
- Compress the target to ignition 
- Inertially confine the target long enough so the fusion yield 
exceeds energy losses (driver efficiency, radiative losses, etc.) 


¢ Equivalent view: Magnetic flux excluded by the conducting liner 
generates external magnetic field pressure, imploding the target 


O-Pinch Concepts 


O-Pinch Basics 


0-pinch coil and circuit 


“Pulsed, high current discharge generates 
a rapidly increasing axial B-field between 
the drive coil and conducting liner 
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peo we * Lorentz force/increasing radial pressure 
Adapted from Miyamoto, K. Plasma Physics for Nuclear 2 . . 
Eason) MIT Press, Combridge, A (1987) implodes a material liner and compresses 


fuel target (typically an FRC plasma) 


- e.g. approach used by Slough et al. in the 
Fusion Driven Rocket (2012-15 NIAC |, II) 


Fusion Driven Rocket 


Driver coils 
“ 


FRC 
vee 


¢ At maximum compression, pressure is 
balanced between the stagnating liner 
. a material, external magnetic field, trapped 
Sy ie Hidon DeNee Hoctecliner Canesten Orter “ADA, internal m agnetic field, and fuel pressure 


2016-4538, 52nd AIAA/SAE/ASEE Joint Propulsion Conf, July 
25-27, 2016, Salt Lake City, UT 
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This Study 


New Approach to Liner Compression 


* Replace the time changing magnetic field generated by the 
pulsed current coil with a target moving rapidly into a steady- 
state magnetic field gradient 


* Equivalent to time changing magnetic field observed in the 
target frame of reference: 
OB 
Vv, — 
2 az 


* Magnetic flux is excluded by/internally trapped by the liner 


[«“B,= 


¢ Unbalanced radial pressure rapidly implodes the fuel target 


¢ Requires implosion time < diffusion time of B-field into liner 
- Typical compression times = several psec, diffusion times = msec 
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Preliminary Concept 


Target Magnetic 
Accelerator Field Coil 
Fuel Target i 
« —*—~ Conversion to Directed Thrust 
ee a ~~_ (magnetic nozzle not shown 
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External field pressure + liner material pressure 
= Internal field pressure + gas kinetic pressure 


Potential Benefits 


How This Helps 


* Replaces pulsed drive coil with steady state superconducting 
magnet, mitigating issues with repetitive, high current pulse 
generation 

- Reduces resistive losses, energy storage requirements 
- Reduces demands on switches, power components, etc. 
- Target can be accelerated to required velocity over longer time period 


¢ Fairly compact linear geometry for in-space applications 
- Strong gradient field produced by small bore magnet 
- Length dictated by target velocity and magnetic field diffusion time 
- Readily incorporates magnetic nozzle for directed plasma thrust 
- Amenable to ground based testing 


* Provides options for axial magnetic field trapping in the target 
and laser preheating to improve target gain 


Challenges and Goals 


Phase | 


¢ Understand the dynamics between the rapidly moving target and 
the gradient field 


- Sufficient pressure generated to provide rapid liner implosion/fusion 
- Growth of liner instabilities (RT, other interchange instabilities; mitigation) 


¢ Initial fuel target design (geometry, liner material, fuel type) 
* Accelerator trades to achieve high target velocities (km’s/s) 


¢ Methods to reduce the required amount of compression 
- Target preheating, embedded magnetic field to reduce thermal transport 


¢ Initial coil design to generate strong magnetic field gradient 
* Evaluate magnetic nozzle to direct the hot expanding plasma 


Pull together into an initial vehicle design and 
comparative mission analysis 


Phase I Study Approach 


Target Fuel 


Reaction Energy Release Cross section vs. center of mass energy 
D+Tsa+n 17.59 MeV 10°; 
} —DT 
D+D—>T+p 4.04 MeV >| —DD 
D+D—3He+n 3.27 MeV 10 |} —p3He 
— pe; 
D+3He+a+p 18.35 MeV oie — a 
a” | ——p'B 
D+°%Li > 2a 2.374MeV | 8 | : 
D+°%Li+p+7Li 5.026 MeV © 10°} 
t 
D+°Lion+7Be 3.38 MeV | 
10 | 
pt+"'B- 3a 8.68 MeV 
[ 
n+%ioTt+a 4.86 MeV 10° = hod = = 
0 1 2 3 4 
n+'Li>T+a+n —2.87MeV 10 10 10 10 10 
E (keV) 


Start with D-T as the “easiest” to ignite for this feasibility study 


Cassibry et al., “Case and Development Path for Fusion Propulsion,” Journal of Spacecraft and Rockets, 52 (2), pp. 595-612, Mar-April 2015 


Approach, cont. 


Accelerator Trades 
Light Gas/Rail Guns 


- Experimentally demonstrated to several km/s 
- Potential erosion of rails/component 
- May be useful for initial ground tests 


¢ Electrothermal Accelerators 
- Investigated to refuel tokomak reactors 
- Ablative plasma arc drives electromagnetic acceleration 
- Experimentally demonstrated to a few km/s with 1-g pellets 


¢ Electromagnetic Accelerators 
- e.g. pulsed inductive macron accelerator = 1-gram at 5-10 km/s 


¢ Laser Acceleration 
- Rapid ablative acceleration of target material 
- Experimentally demonstrated to several 10’s of km/s 


Approach, cont. 


Reduce Compression Requirements 


¢ Fuel preheating reduces amount of adiabatic compression 
required to reach fusion temperature conditions: 


n\ 2 eas 
T.=T ‘0 rg = initial fuel radius 
f e ry r= final fuel radius 
- Option: Laser provides both ablative target acceleration and fuel 


heating prior to compression in magnetic field gradient 


¢ Initial magnetic flux trapped within fuel provides strong axial 
field when compressed: 


n 2 
By =B (2) 
f 
oO ry 


- Reduces radial ion and electron thermal transport losses 
- Enhances compressed fuel heating by a-particle energy deposition 


ee 


Approach, cont. 


NASA Inno 


Analytic Models 


¢ MATLAB modeling for preliminary field calculations, adiabatic liner 
compression, target sizing estimates; using to guide initial inputs 
to comprehensive 3D numerical model (SPFMax) 


¢ Adapting semi-analytic MIF model* to assist with rapid trade 
space evaluations and target performance estimates 

- Liner dynamics and compression, fuel ionization and adiabatic heating, 
optional fuel preheating (laser) absorption), a-particle energy deposition, 
radiative losses, ion and electron thermal conduction losses, end losses, 
magnetic flux losses, fusion reaction rates, energy balance and gain 
calculations 

- Reasonable published comparisons with LASNEX numerical models 

- Modifying equations to incorporate high speed target interactions with 
gradient magnetic field and corresponding liner compression dynamics 


*McBride, R. and Slutz, S., “A semi-analytic model of magnetized liner inertial fusion,” Physics of Plasmas 22, 052708 (2015) 


Approach, cont. 


Numerical Model: SPFMax 


Smooth Particle Hydrodynamics with Maxwell equation solver* 


* 3D multiphysics code, vectorized and GPU-enabled in Matlab; 
developed by UAH for MIF and fission/fusion hybrid research 


* Includes tabular equations of state to model variable levels of 
ionization, thermal conduction, radiation emission and 
absorption, shock capturing, real viscosity, electromagnetic field 
propagation and forces in the plasma, a self-consistent circuit 
model, and nonlocal absorption of fusion ion product energy 


* Resolves vacuum/plasma interface and nonlocal transport of 
charged particle deposition 


* Currently completing the integration and testing of the SPFMax 


electromagnetic equation solver 
*Additional details in backup slides 


Approach, cont. 


Preliminary SPFMax 
Results 
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Initial Magnetic Field Topology in 
Converging/Diverging Nozzle 
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Density slice and surface plot during initial motion of target through 
entrance of coil. Azimuthal current induced on target by vxB back EMF 


Approach, cont. 


Mission Modeling 


¢ Initial parametric assessments using approximation techniques 
to optimize performance requirements 
- e.g. Gilland, J. H., “Mission and System Optimization of Nuclear Electric 
Propulsion Vehicles for Lunar and Mars Missions,” NASA CR 189058 
* Model inputs: power, specific impulse (Isp), thrust or efficiency 
vs. Isp, system mass, desired payload mass, propellant 


* Initial cases being evaluated include round-trip Earth-Mars and 
round-trip Earth-Saturn 


* Detailed mission analysis using Copernicus code once engine 
performance and vehicle design are more refined 


- Copernicus Trajectory Design and Optimization System software 
(https://software.nasa.gov/software/MSC-25863-1; available for use by federal 
employees and contractors to the federal government ; installed at MSFC and OAI) 


Initial Mission Modeling Results 
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Test case: Payload mass delivered 
to Mars with NIAC PUFF concept 
(2013/Adams et al.) 

* 16 kW/kg, 1 GW pulsed system, 
evaluated over a range of Isp from 
10,000 s to 50,000 s 

¢ Baseline the approximation method 

¢ Adapt for gradient-field concept 
performance parameters 


Gradient field concept: Preliminary modeling estimates indicate 

one-way trip times of less than 1-year for both Mars and Saturn: 
* Mars: 60-90 days, 100-150 mT payload from LEO, 2 SLS-Block 2 launches 
¢ Saturn: 200 days, 50-100 mT payload from LEO, 2 SLS-Block 2 launches* 
¢ Initial results indicate Isp > 20,000 s provides greatest benefit 


*Jackman, A., “Space Launch System” (2016) https://sites.nationalacademies.org/cs/groups/ssbsite/documents/webpage/ssb_171779.pdf 
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Expanding Plasma 


Approach, cont. 
Initial Vehicle Concept 
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Ignited Plasma 


Mass Breakout 
Dry Mass 94 mT 


(30%) MGL 40 mT 
Fuel 56 mT 
IMLEO 240 mT 


S=s 
SS 
SSS 
SSS 
SS 
ss 
Fuel Tanks 
Initial breakout based on previous 
Fusion concepts at MSFC 


Deep Space Habitat 


Summary 


* Developing the tools needed to evaluate concept performance 
in a mission context: 
- Analytic models to provide initial performance estimates 
- SPFMax code for detailed understanding of target implosion physics, 
energy gain, and directed plasma exhaust 
- Trajectory approximation methods to bound mission parameters 
- NASA Copernicus trajectory code for more detailed mission analysis 
- Updating prior MSFC and UAH fusion vehicle analysis with new engine 
design and performance parameters 
* Incorporates steady advancements in larger DOE terrestrial 
magnetoinertial fusion experiments 
- Guides the target designs, implosion conditions, and magnetic field 
strengths that need to be reproduced in the gradient-field concept 
- Innovative MIF approach adapted for in-space propulsion applications 
- Allows continued progress toward mission-enabling fusion vehicle design 
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Primer on Smooth Particle Hydrodynamic Method 


* Gridless Lagrangian technique wit 


approximations 
1. Integral approximation 


A,(r)={ A(r')W(r—r',h)dr' 


2. Integral approximated by 
summation interpolant: 


A= A,V,We(t—r',h) 
b 


Particle and 


* Derivatives are obtained by: 


VA,= >, A,V,V Walt—r'sh) 
b 


60 nearest 


sigleievo neighbors 


: i000 particles 


x (m) 
60 jets in 
a soccer-ball 
configuration 


z(m) 
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Approach motivated by 


Fast turnaround on 
3D simulations 
¢ Mixing of particles 


¢ 0 
a . 


y(m) x (m) 


SPFMax, cont. 


Governing plasma equations in SPFMax: 


Ou 1 

Bp pV PtV +H 4XB 

de, P; T 

apap utp Vu-V(kVT))+0, 
Oe 
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SPFMax, cont. 
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Implementation of the Electromagnetic solver uses an external lumped 
element circuit model, transmission line network for the SPH particles, 
and scalar and magnetic vector potential for the electromagnetic fields 
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SPFMax, cont. 


SPFMax Radiation Matter Interactions 


Magneto-inertial Fusion Target 
¢ Radiation/matter interaction using Rays 


DT Fuel * Nonlocal deposition of multigroup 
Target radiation using absorption opacities 


* Nonlocal deposition of charged particle 
energy using stopping power 


Radiation/matter properties, like 
opacities, are interpolated to ray 
segments based on local particle values 


0.05 , 


Argon 
Liner . 


Rays Segment 


a ie 
Radiation 
source 
* Power is attenuated as it propagates 

through the branch 


* Power is absorbed in particles based on 
Rays Branches local ray segment property, particlessize, and 
distance from radiation source 


GRADIENT FIELD IMPLODING LINER FUSION PROPULSION SYSTEM 
A NEW APPROACH TO MAGNETOINERTIAL FUSION FOR IN-SPACE PROPULSION 
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BACKGROUND: Magnetoinertial fusion (MIF) employs a high, rapidly pulsed current in a target liner Why Fusion? 
to generate an intense magnetic field, creating an unbalanced pressure that radially implodes the socket Patton eoolapevstenaetnaioner cau 
target to achieve fusion conditions. Compressed magnetic fields internal to the fuel target reduce m r Av = 104— 10° m/s 
. 7 < 7 . . . — Vv e i 
electron thermal diffusion losses and confine fusion generated a-particles to improve fuel heating. ut —e /|ve Reg onabie PaylOadinass actors 
require propellant exhaust velocities 
Mo comparable to mission Av 
INNOVATION: We replace the pulsed current coil and stationary target used in MIF O-pinch concepts m,=initial mass Av=velocity change _—* Fusion offers high specific power 
m,=finalmass __v, = exhaust velocity (kW/kg) to reduce trip times 


with a fast moving target fired axially into a static, high gradient magnetic field. This is equivalent to 
generating an intense time changing magnetic field in the target frame of reference, resulting in an 


unbalanced radial pressure that rapidly implodes the fuel target. MIF Lind!-Widner Diagram (DT cylinder at stagnation)* 


Advantages of MIF 
¢ Adiabatic radial compression 
generated by unbalanced 
magnetic field external to target 
* Compresses a magnetic field 


Q-pinch coil and circuit for Compression Steady-State Gradient Magnetic Field for Compression 


Target 


Coil 


Target 
Accelerator 


Magnetic 
Field Coil 


Fuel Target 


Ss 
oO 
| ——> ra — Conversion to Directed Thrust J : 
aici tas_aaa cael =a a SIuiegnetie havin bor towh) © internal to the target, reduces 
= electron thermal losses and 
Vacuum = enhances alpha particle self- 
Initial target wall P " 2 heating of the fusion fuel 
diameter : \ = 


Energy storage and 
pulse forming network 


Adapted from Miyamoto, K. Plasma Physics for Nuclear 
Fusion, MIT Press, Cambridge, MA (1987) 


¢ Reduces the areal density (pr) 
threshold for ignition to well 
below conventional Inertial 
Confinement Fusion (ICF) values 


Compress Burn Expand 


1 
———_— 0.001 0.01 0.1 1 
. : : *Basko, M. et al., “Ignition conditions for magnetized 
. Pulsed, high current discharge generates a rapidly . . ; fuel pR (g/cm’) target fusion in cylindrical geometry,” Nuclear Fusion, 
é - : : : ‘ External field pressure + liner material pressure 
Increasing axial B-field between the d rive coil a nd = Internal field pressure + gas kinetic pressure Vol. 40, No. 1 (2000) 
conducting liner 5d OB 
: : : ; 1x B, = vy 
¢ Increasing radial pressure due to time changing dz ise _— — 
magnetic field implodes a material liner and ¢ Magnetic flux is excluded by/internally trapped by the liner | 
compresses the fuel target ¢ Unbalanced radial pressure rapidly implodes the fuel target a _ ne 


0.05 


BENEFITS: 
e Replaces the pulsed current coil and PFN with a steady state superconducting magnet 
- Mitigates associated issues with rapid, repetitive, high current pulse generation in the drive coil 
- Target acceleration over a longer time period reduces demands on switches, power components, etc. 


¢ Preliminary analysis indicates a fairly compact engine geometry 
- Strong gradient magnetic field is produced by a constant current, small bore coil 


- Maximum system length is dictated by target velocity and magnetic field diffusion time yim) S08 Ly ss SS ig Oe 
¢ Linear geometry provides options for generating/trapping an initial axial magnetic field in the SPFMax: Density slice and surface plot during initial motion of target through 
target and for pulsed laser preheating to improve fuel ignition and compression temperatures entrance of coil. Azimuthal current is induced on target by vxB back EMF 
KEY PHASE | GOALS: Preliminary Vehicle Design Preliminary Mars Mission Analysis 
e Establish requirements for magnetic field coil, target accelerator, and Magnetic Nozale Expanding Plasma Isp (log) 10 fedtcad 
encapsulated fuel target, leading to an initial system concept design LS (mT) 


¢ Perform analysis and detailed SPFMax modeling to understand the - Teo 
dynamics between the rapidly moving target and the gradient field ee 100 
- Sufficient radial force to provide rapid liner implosion and fuel ignition Mass Breakout il | % 


(30%) MGL 40 mT 200 
Fuel 56 mT 
IMLEO 240 mT 


- Growth of liner instabilities (interchange instabilities and mitigation) 


¢ Develop integrated vehicle concept based on prior MSFC/UAH designs 


100 |. 
50 


100 
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¢ Model representative Earth-Mars and Earth-Saturn round-trip missions 


150 = One-way Earth to Mars trip 
for 201s NIAC PUFF concept; 


Taliat-] Mm o)g=t-] colli mel-ky-1e Mela ol g-V (e013 


using NASA Copernicus trajectory code to establish concept benefits FuslelN eee SUE (days) equivalent length method 
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